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1.1 The Thesis Question 

Can air and ground resonance stability be improved by judicious use of composite 
couplings in the wing and kinematic couplings in the blades of a soft-inplane hingeless tiltrotor? 

1.2 Background and Motivation 

The tiltrotor offers the best mix of hovering and cruise flight of any of the current V/STOL 
configurations. One possible improvement on the tiltrotors of today designs would be using a 
soft-inplane hingeless hub. The advantages to a soft-inplane hingeless hub range from reduced 
weight and maintenance to reduced vibration and loads. However, soft-inplane rotor systems are 
inherently in danger of the aeromechanical instabilities of ground and air resonance. Furthermore 
tiltrotors can be subject to whirl flutter. At least in part because of the potential for air and ground 
resonance in a soft-inplane rotor, the Bell XV-15, the Bell-Boeing V-22 Osprey, and the new Bell 
Augusta 609 have stiff-inplane, gimballed rotors which do not experience these instabilities. In 
order to design soft-inplane V/STOL aircraft that do not experience ground or air resonance, it is 
important to be able to predict these instabilities accurately. Much of the research studying the 
stability of tiltrotors has been focused on the understanding and prediction of whirl flutter. As 
this instability is increasingly well understood, air and ground resonance for a tiltrotor need to be 
investigated. 

Once we understand the problems of air and ground resonance in a tiltrotor, we must look 
for solutions to these instabilities. Other researchers have found composite or kinematic couplings 
in the blades of a helicopter helpful for ground and air resonance stability. Tiltrotor research has 
shown composite couplings in the wing to be helpful for whirl flutter. Therefore, this project will 
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undertake to model ground and air resonance of a soft-inplane hingeless tiltrotor to understand 
the mechanisms involved and to evaluate whether aeroelastic couplings in the wing or kinematic 
couplings in the blades would aid in stabilizing these instabilities in a tiltrotor 

A brief history of tiltrotors and the literature related to this project is given in the next 
section. The remaining sections provide detail about the current analysis. The modeling techniques 
used to model the blades, hub, rigid pylon, and elastic wing are described in detail. The validation 
study used to give confidence to the model is shown with comparisons to a soft-inplane rotor and 
to a gimballed rotor Some results from the parametric studies to look at the couplings in the blades 
are given followed by a discussion of the future research needed and the time line for completing 
that research. 



2 Literature Review 
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2.1 Tiltrotor History 

One of the biggest hurdles facing a tiltrotor designer today results from the relatively 
few number of tiltrotors that have been built: the knowledge that comes from experience is very 
limited. It behooves us to spend a little time remembering what has been done before embarking 
anew. 

The first tiltrotor built and flown was the Transcendental Model 1-G as shown in Fig. 1 
[16] . The Model 1-G was developed without governmental support of any kind, a rarity today, and 
was the first demonstration of the capabilities of a tiltrotor configuration. The Model 1-G achieved 
first free flight on July 6, 1954 at Bellanca Field, New Castle, Delaware, and first conversion 
to airplane mode in December, 1954. The Model 1-G used three-bladed, fully-articulated rotors 
(i.e., rotors with flap and lag hinges and pitch bearings-see Fig. 2). After 23 hours of air time 
over more than 100 flights, the Model 1-G crashed on July 20, 1955 and was not rebuilt [28] . In 
1956 the Air Force awarded Transcendental a contract to build another tiltrotor. The result was the 
Transcendental Model 2-G. basically the same design as the Model 1-G but with more powerful 
engines. The Model 2-G was never flown [35] . 

Bell Helicopters then developed the Bell XV-3, one of two designs chosen from a Air 
Force / Army “convertiplane” competition in 1951 where 17 companies submitted 19 designs [28] 
[35] . The XV-3 also had three blades and a fully-articulated hub. First flight for the XV-3, the 
second tiltrotor to fly successfully, was in 1955. The first XV-3 built crashed, and a second XV-3 
was built but this time with two-bladed teetering rotors (i.e., rotors with two blades rigidly attached 
to one another with the shaft as a pivot point for flapping-see Fig. 3). The XV-3 led researchers 
to discover and solve many technical challenges of the tiltrotor configuration [16] . In 1968, Bell 
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began to develop the Model 300, now three-bladed but with gimballed rotors (i.e., rotors with 
three blades rigidly attached to a hub which is free to pivot about the shaft on a universal joint). 
Dynamic rotor/pylon instability tests were conducted, but the vehicle was never built. The stability 
margins which had been a problem with the articulated rotors did not seem to be a problem with 
this gimballed design: the Model 300 was stable well beyond its design dive speed [35] . 

NASA and the U.S. Army initiated a proof-of-concept program in April, 1973. The 
contract was awarded to Bell Helicopter to construct what began life as the Model 301 but became 
the XV-15 Tiltrotor Research Aircraft shown in Fig. 4. The XV-15 has the best safety record of 
any tiltrotor built to date. The 6.5 ton aircraft flew up to 346 mph at 16,000 ft. in June, 1980. It 
has been called the “most successful high speed V/STOL ever built” [16] . 

In December 1981, with the XV-15 tests progressing apace, the Deputy Secretary of 
Defense initiated a development program for a single aircraft which would serve the U.S. Marine 
Corps, the U.S. Air Force, and the U.S. Army. Budgets were shrinking so a single aircraft, or 
at least several versions of a single aircraft, had the potential to save the Defense Department 
significant sums of money. The requirements were stringent: for example, the craft must hover 
at high altitudes at high temperatures with 8300 lb. payload. It must also be able to self-deploy 
anywhere in the world without in-flight refueling [27] . Conventional helicopters, compound 
helicopters, tiltrotors, and lift fan concepts were considered to meet these goals, but the tiltrotor 
was chosen as the best configuration to meet the various requirements. A contract was awarded 
to the team of Bell Helicopter Textron and Boeing Helicopters in 1982 with a view to develop 
and produce the V-22 Osprey (Fig. 5), originally the JVX, eventually to deliver almost 1000 
tiltrotors for four branches of the armed forces [27] . The first flight occurred on March 19, 1989. 
Development of the V-22 is continuing at this time: the second wave of test vehicles has flown 


approximately 1400 hours to date. 
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Several other tiltrotors have reached the design or testing stage, including the Bell 266, 
the Grumman Helicat, and the Boeing 222. Each offered valuable information about the tiltrotor 
configuration, but none of these aircraft ever flew [22] , [26] . 

2.2 Advantages of Tiltrotors 

The primary advantage of tiltrotors over helicopters is speed. Because the wing offloads 
the rotor in forward flight, the rotor can be used to provide forward speed without also providing 
lift. It is estimated that a tiltrotor can reach a rescue site on average 2-3 times faster than present 
rescue aircraft [35] . Also, the limitation of edgewise flight which a helicopter faces is avoided by 
tilting the rotor, in effect placing the rotor into a situation akin to climb in a helicopter [35] . But 
if the only desire is speed, fixed wing aircraft are much superior. The advantage of tiltrotors over 
fixed wing is obvious: they can hover, albeit not as efficiently as helicopters. The advantage of 
tiltrotors over a Harrier or other jet V/STOL is lower disk loadings, or the ratio of thrust to disk 
area (Fig. 6). This lower disk loading makes the tiltrotor much more efficient in hover (since 
induced power loss is proportional to thrust multiplied by downwash velocity) and affords a lower 
downwash velocity and lower tip-speed than other alternatives. As the tip-speed is lowered, the 
noise level is decreased. As the downwash velocity is lowered, the tiltrotor is less bothersome to 
people on the ground. The tiltrotor is an ideal choice when both speed and hovering capabilities are 
important, though clearly it will never replace the helicopter for a truly efficient hovering machine 
or the fixed wing for long fast flights from prepared runways. 

There are other possibilities for the mix between speed and hovering ability, many of which 
can be considered compound helicopters. But a study performed by the Group for Aeronautical 
Research and Technology in Europe (GARTEUR) undertaken in the late 1980’s compared several 
compound helicopters with a tiltrotor and concluded that the tiltrotor was more economically 



feasible for a large range of operations than the compound helicopter [11] . 
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2.3 Soft-/Stiff-Inplane Hingeless Tiltrotors 

Tiltrotors with stiff-inplane, gimballed rotors have shown their potential in the XV- 15 
which is still flying, the V-22 which is almost through with flight testing, and the Bell Augusta 
609 which is available for ordering. Clearly this rotor system is effective. Why then do we want to 
reconsider the design? 

The tiltrotor offers intriguing possibilities for civil transport. For example, the “Civil 
Tiltrotor Missions and Applications” research study concluded that tiltrotors could capture a 
major share of the travel in the northeast corridor provided vertiports were available [48] . 
NASAs Civil Tiltrotor Initiative is dedicated to helping this become a reality if possible. To be 
economically viable, a civil tiltrotor needs to be as inexpensive and efficient as possible. One 
possible improvement would be to use a soft-inplane hingeless hub. Soft-inplane rotor systems 
have an inplane natural frequency lower than the rotor speed. This is the norm for articulated rotors 
where a typical lag frequency is around .2 or .3/rev. But as hingeless rotors are increasingly chosen 
for new designs, the V/STOL designer has a choice: the lag frequency shouldn't be too close to 
1/rev because of excessive blade loads but it can be either above or below 1/rev. Helicopters have 
found advantages in using soft-inplane hingeless rotor systems such as reduced vibration and loads 
and potentially improved handling qualities. There are fewer parts in a hingeless hub resulting 
in potentially less weight and improved reliability. The Comanche is a prime example of these 
advantages [38] . 

But it remains true that the most successful tiltrotors, the XV- 15, the V-22, and the new 
609, all have stiff-inplane, gimballed rotors. This is at least in part because a soft-inplane rotor is 
inherently in danger of the instabilities of ground and air resonance (see [18] ). Ground resonance 
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is a catastrophic instability which can destroy aircraft and anything inside them. Air resonance is a 
less catastrophic instability which causes large vibrations and severely degraded handling qualities. 
The challenge of building an aircraft which satisfies all the mission requirements and doesn’t suffer 
from air and ground resonance is more complicated for tiltrotors than it is for helicopters since 
tiltrotors also suffer from whirl flutter [20] . 

In order to design soft-inplane V/STOL aircraft that do not experience ground or air 
resonance, it is important to be able to predict these instabilities accurately. Analytical computer 
modeling to predict air and ground resonance in soft-inplane helicopters has been available for 
some time. Much of the research studying the stability of tiltrotors has aimed to understand and 
predict whirl flutter. As this instability is increasingly well understood, air and ground resonance 
for a tiltrotor need to be investigated. 

Once these instabilities in a tiltrotor are well understood, the question becomes what can 
be done to ensure stability. Many studies have found composite or kinematic couplings in the 
blades of a helicopter helpful for ground and air resonance. Tiltrotor research has shown composite 
couplings in the wing to be helpful for whirl flutter. Therefore, this project will undertake to model 
ground and air resonance of a soft-inplane hingeless tiltrotor to better understand the mechanisms 
involved and to evaluate whether aeroelastic couplings in the wing or kinematic couplings in the 
blades would aid in stabilizing these instabilities in a tiltrotor. 

2.4 Instabilities 

One of the great challenges in aircraft design is to be very sure that there are no instabilities 
encountered within the operating regime of the aircraft. The tiltrotor has several possible 
instabilities which must be avoided. Whirl flutter is a concern for all tiltrotors while air and ground 
resonance concern only soft-in-plane tiltrotors. 



2.4.1 Air and Ground Resonance in a Helicopter 
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Ground resonance was first observed in the autogiros in the 1930’s and early 1940’s. This 
instability was not well understood at that time and was therefore difficult to avoid by design. 
The American Helicopter Society listed the characterization of air and ground resonance as one 
of the top twenty advances in vertical flight history [7] . Ground resonance occurs when the low 
frequency lag mode of the rotor (in the non-rotating frame) coalesces with a fuselage or body 
frequency. The cyclic lag modes £ lc and ( ls produce a wobble of the rotor center of gravity. This 
shifting fore and aft and laterally of the effective mass of the rotor can couple with the vibration 
of the fuselage to produce huge displacements. On a helicopter sitting on the ground, the body 
frequencies are governed primarily by the stiffness of the landing gear-aerodynamic forces play 
very little role. Ground resonance can destroy an aircraft in seconds (as shown by the Kellet XR-2 
shown in Fig. 7). 

Air resonance is a similar aeromechanical instability mathematically but occurs in the air, 
often in a hingeless rotor [18] . Unlike for ground resonance, aerodynamic forces play an important 
role in air resonance. Air resonance is rarely catastrophic, but since this instability causes very 
large vibrations (which leads to issues of passenger comfort and material fatigue life) and much 
poorer handing qualities, regions of instability should be avoided if possible [37] . 

The instability of ground resonance was first explained mathematically by Coleman and 
Wagner independently in 1943 [9] . These early analyses neglected aerodynamics, certainly 
a good first approximation for ground resonance. Both Coleman and Wagner developed a 
simple four-degree-of-freedom model which included high and low frequency lag modes (in the 
non-rotating frame) for the blade and body pitch and roll. The graph shown in Fig. 8 shows the 
rotor rotation speed on the horizontal axis and the frequencies of the various degrees of freedom 
on the vertical axis. As the low frequency lag mode coalesces with the body mode an instability 



occurs. This form of graph has come to be called a Coleman diagram [9] . 
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The damping associated with this simple model is quite distinctive as shown in Fig. 9. 
This plot shows the damping ratios from a four-degree-of-freedom system, blade cyclic lag plus 
body lateral and longitudinal motion. With no aerodynamics present, the damping of all four 
degrees of freedom is zero except in the regions of instability. In these regions the low frequency 
lag mode has coalesced with one of the body frequencies: one mode becomes unstable and the 
other mode becomes even more stable. In a helicopter it is the low frequency lag mode which goes 
unstable. As aerodynamic forces are added, the damping is no longer zero outside of the unstable 
regions. When this is the case one of the modes will be more damped than the other. It is typically 
the less damped mode that becomes unstable. Again, in a helicopter this usually means that the low 
frequency lag mode becomes unstable while the fuselage mode becomes even more stable since 
the rotor body modes are usually more damped than the rotor lag modes. Preliminary results with 
this tiltrotor analysis show that it is typically the wing modes that go unstable. 

Coleman's formulation gives rise to a very simple stability criterion shown below: a 
helicopter would be stable if 


r * r , (1 f“y\ 2 

C < ( " > 


\nJ m* 


where is is half the ratio of the rotor mass to the support mass, ^ is the non-dimensional body 
lateral mode, C£ is the damping in the lag mode, C* is the damping in the body lateral mode, and 
is the rotor lag frequency (assuming that the frequency for the body mode is not the same for 
both lateral and longitudinal motion) [8] . This simple criterion shows that it is necessary to have 
damping in both the body mode and in the lag mode. Thus landing gear dampers and rotor lag 
dampers exist on almost all helicopters flying today. The presence of the ratio of the rotor mass 
to the support mass is instructive: if the rotor mass is large with respect to the body mass, more 
damping is needed to stabilize the vehicle. Also, this simple stability criterion makes it clear that a 



12 

stiff-inplane rotor will not suffer ground resonance since (1 — uq ) is negative [18] . 

Coleman's basic theory was not enough to eliminate the problems with other helicopters, 
including the Sikorsky XR-4, R-5, and R-6 [7] . Clearly a more exact analysis of this phenomena 
is desired since avoiding ground resonance is a prime design requirement. Advances in ground 
resonance prediction have come along steadily since the 1970’s with a plethora of studies. In 1970 
Lytwyn, Miao, and Woitsch expanded the basic analysis to hingeless rotors where air resonance is 
more of a problem [25] . Dynamic inflow was shown to be important in the early 1980’s starting 
with Johnson in 1982 and Friedmann and \fenkatesan in 1986 [19] , [12] . The importance of an 
accurate trim model to air resonance stability was shown by Nagabhushanam and Gaonkar and 
later by others [31] , [30] . Ground and air resonance stability is now incorporated into many 
comprehensive rotorcraft codes. Some additional contributions will be discussed in the next 
section. 

Comprehensive codes to predict ground and air resonance in tiltrotors have been slower in 
coming. The lack of reliable methods of ensuring ground resonance for a tiltrotor has certainly 
contributed to the pervasive use of gimbal technology in the tiltrotors flying today Going to a 
stiff-inplane gimballed hub for the XV- 15 and the V-22 eliminated these problems. 

2.4.2 Whirl Flutter in a Tiltrotor 

In 1960, two Lockheed Electra aircraft crashed with loss of life in what is now called 
propeller/nacelle whirl flutter or simply whirl flutter [20] . Whirl flutter is a dynamic instability 
where the inplane forces couple with wing motions to become unstable. This possibility was 
first published by Taylor and Browne in 1938 and later investigated academically by Scanlan 
in the 1950’s [47] , [43] . The cases of whirl flutter in the Electra, occurring as they did in a 
propeller-driven, fixed-wing aircraft, suggested the possibility that the tiltrotor configuration might 
also be subject to a similar instability [42] . 
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In a 1962 wind-tunnel test of the XV-3 in the NASA Ames full-scale tunnel, an instability 
driven by a coupling between the proprotor and the pylon was discovered [21] . This finding led 
many researchers to look into whirl flutter in tiltrotors. Kvatemik reviewed the state of the art of 
whirl flutter research in 1976 and again in 1991 and 1992 [22] , [23] , [24] . 

Whereas a typical ground resonance plot is with a rotor speed sweep, a typical whirl flutter 
plot is versus forward speed as is shown in Fig. 10. The wing modes marked c, b , and t in the plots 
here (chordwise bending, beamwise bending, and torsion) all go unstable at high enough forward 
speeds. This instability is very violent and can destroy the aircraft. 

2.4.3 Air and Ground Resonance in a Tiltrotor 

Twenty-six years ago Boeing conducted two wind tunnel tests for the Model 222 which 
was designed as soft-inplane hingeless tiltrotor. The first test in 1973 used a 1/4.6 Froude-scaled 
Model 222 in the Boeing V/STOL wind tunnel. Steady state rotor and aircraft loads, performance 
data, and stability derivatives from mid-transition to cruise were obtained [29] . The second 
test was a full-size, semi-span model test in the NASA Ames 40’ by 80’ wind tunnel. This test 
identified regions with degraded stability and investigated some active feedback systems to add 
damping to the lightly damped wing modes. All the tests done in the wind tunnel were in cruise 
mode; no ground resonance data was gathered [26] . This test was compared to the 1971 Bell 
Helicopter test of a 25-foot gimballed rotor in the same NASA 40’ by 80’ wind tunnel [17] . Both 
of these tests used the same NASA wing and pylon making it easier to compare them. 

2.5 Possible Solutions 

During the years from 1970 to 1998, advances in computing power developed to allow 
more exact characterizations of the instabilities of ground resonance, air resonance, and whirl 
flutter. The obvious mechanisms of ensuring the stability of the V/STOL vehicle were to use 
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either stiff-inplane rotor systems or lag and landing gear dampers. These solutions neglected the 
advantages of a soft-inplane hub and added weight and complexity. The search was then for other 
ways to stabilize these instabilities. 

One possible solution to the stability problems of air resonance, ground resonance, and 
whirl flutter, is the use of aeroelastic couplings. Composite tailoring techniques can now introduce 
couplings into wings or blades. For example, consider a composite beam with two layers, one with 
fibers running at 45 degrees and one with fibers running down the length of the beam. This kind of 
unbalanced laminate structure will exhibit bending-torsion coupling. Pushing up on the end of the 
beam will twist the beam. This is the sort of construction that could be used in the spar of a tiltrotor 
wing or blade (Fig. 11) or in the flexbeam of a rotor blade (Fig. 12). 

Kinematic couplings in the blades could be created in several ways other than with 
composite tailoring. Skewed hinges in an articulated blade can create pitch-lag or pitch-flap 
couplings (Fig. 13). The order of the hinges when they aren’t coincident create couplings from 
the steady coning or lag angles. Feedback systems can couple pitch and blade motion. Also, a 
pitch-link outboard of a flexure as in this example causes pitch-flap couplings (Fig. 14). 

2.5.1 Blade Couplings For Helicopter Air and Ground Resonance 

Using kinematic couplings in the blades to stabilize air and ground resonance for a 
helicopter is not an uncommon concept in the literature. Burkam and Miao started this research in 
1972 [ 6 ] . In 1977 Ormiston published a landmark paper investigating many different contributing 
factors. He discovered that combinations of pitch-flap and pitch-lag couplings with structural 
flap-lag coupling were most helpful for stabilizing air and ground resonance. Ormiston also 
showed that ground resonance was very different from air resonance and that the factors which 
improved air resonance stability margins (for example pitch-lag coupling, £4 < 0 , where the blade 
pitches nose-down as it lags back) destabilized ground resonance [36] . 
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Bousman worked with these couplings experimentally in 1981; this research showed that 
pitch-lag coupling (64 < 0 ) improved stability away from the critical regions of air resonance 
instability but was less helpful immediately around the instability [5] . In 1993 Smith and Chopra 
investigated composite tailoring of the blades to postpone or eliminate ground and air resonance 
instabilities for a helicopter [45] . Composite tailored blades with pitch-lag coupling improved air 
resonance stability and pitch-flap couplings reduced the air resonance stability margins though 
this effect was less powerful than the pitch-lag couplings. Gandhi and Hathaway published an 
optimization study that showed that pitch-lag with £4 < 0 , pitch-flap coupling where the blade 
pitches nose-down as it flaps up (63 > 0), and a lag-pitch-flap hinge sequence formed a good 
combination to stabilize ground resonance [14] . 

2.5.2 Blade Couplings For Tiltrotor Whirl Flutter 

Gaffey published a study in 1969 which considered pitch-flap coupling for whirl flutter 
in a tiltrotor. Gaffey found that pitch-flap coupling with 6 3 < Ohad a stabilizing trend on whirl 
flutter and that 63 > 0 destabilized whirl flutter [13] . More than twenty years later, Nixon 
showed that elastic extension-twist coupling of the rotor blade can improve tiltrotor aerodynamic 
performance but that this coupling degrades the stability for whirl flutter [35] . Nixon's research 
also investigated composite couplings in the blades to stabilize whirl flutter. This study found 
that elastic bending-twist coupling in the blades has potential to increase whirl flutter stability 
and that the amount of coupling required is feasible for performance and blade loads. Fig. 15 
shows an increase in whirl flutter stability based on flatwise bending-twist couplings in the blade 
taken from Ref. [35] . Here A = is a measure of the off-diagonal coupling between the 

flatwise-bending and twist of the blade. A is sensitive to the forward velocity unlike kinematic 
pitch-lag coupling because as V increases the blade has more collective pitch so that more of 
the flatwise bending in the lag direction. As V increases from 100 to 400 knots, a A value of .1 



increases from a 64 value of 0.4 to 0.7. 
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2.5.3 Composite Tailoring of the Wing for Whirl Flutter 

In 1995 a group of researchers from Bell published an analysis showing that tailoring 
the wing with composites could allow the V-22 to use an 18% thickness-to-chord ratio wing 
instead of a 23% ratio without degrading the whirl flutter stability margin [39] . Using comparable 
thicknesses, a 24 knot increase in flutter speed was predicted. In 1997 a wind tunnel test at NASA 
Langley showed a 30 knot delay in whirl flutter onset (Fig. 16) [10] . The 18% thick wing used 
wing beamwise bending/torsion coupling in the skin to decrease the natural coupling between 
bending and torsion due to the offset of the pylon and rotor masses from the pylon attachment 
point. This improved the stability of the wing bending mode but degraded the stability of the wing 
chord mode. Tailoring the wing spar caps alleviated the loss of stability for the wing chord mode. 

Barkai and Rand showed that chordwise bending-twist in the wing was beneficial for whirl 
flutter stability in gimballed, stiff-inplane tiltrotors [1] . Fig. 17 shows the increase in whirl flutter 
stability she calculated using wing chordwise-bending-twist coupling. In this figure 2m — 2 N e~ N 
where Ng is the number of plies oriented at an angle 9 and N is the total number of plies; C\q 
is the off-diagonal coupling modulus for bending-twist coupling in the composite layup [41] . 
Though this study and the Bell study discussed above were for gimballed, stiff-inplane tiltrotors, 
their importance lies in a stability increase for a tiltrotor using aeroelastic couplings in the wing. 

2.6 Modeling 

Aerospace engineering has come a long way from the days when the an aircraft could be 
built and flown in a year. Stringent requirements are placed on new aircraft today for performance, 
efficiency, and safety. Nor are funds so free-flowing that we can afford to tread far down the road 
for one design only to find it untenable. Comprehensive modeling codes are used to test as much as 
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possible for each potential design as inexpensively as possible. Accurate modeling of the tiltrotor 
is essential to predicting its stability. Since the XV-3, the XV-15, and the V-22 are all stiff-inplane 
gimballed rotors, the vast majority of studies conducted for tiltrotor stability have focused on whirl 
flutter analysis. 

Aeroelastic analysis began in the 1960’s and expanded significantly in the 1970’s. In 1974 
Johnson compared the Bell and Boeing rotors discussed above using an analysis that included 
hover, conversion, and forward flight and later expanded this model to a comprehensive aeroelastic 
model which became CAMRAD [17] . In 1984 and 1985 Kvatemik developed the Proprotor 
Aeroelastic Stability Analysis (FASTA) code at NASA Langley [24] . By the third version of 
PASTA eleven degrees of freedom were modeled including rigid flapping and lead-lag motion of 
the blade, gimbal motion, and a rigid control system. This code as many others of that era included 
only cruise mode. Also, a comprehensive aeroelastic analysis called DYN4 (for rigid blades or 
DYN5 for elastic blades) has been developed by Bell Helicopters [46] . 

In 1993 the University of Maryland Advanced Rotorcraft Code (UMARC) was expanded 
to allow whirl flutter research. This analysis, which included hover, transition, and forward flight 
modes, consisted of a gimballed hub with blades and wing modeled using the finite element 
method [35] , [3] . 

Srinivas expanded the UMARC code in 1995 to include the entire tiltrotor including the 
whole wing span and both rotors using the finite element method (Fig. 18) [46] .One of Srinivas’ 
conclusions was that for whirl flutter the symmetric modes modeled using a wind tunnel, semi-span 
model were the modes that became unstable first and hence set the flutter speed. This indicates 
that a good first step in modeling a tiltrotor would be to use a semi-span model. 

In 1995, a group of researchers from the Technion Institute in Israel presented a new 
analysis based upon a numerical technique to preserve the symbolic exactness of the equations of 
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motion. This method retained the partial derivatives for the stability analysis and allowed them to 
discuss some of the causes for whirl flutter in tiltrotors [40] ,[41]. 

2.7 Conclusions 

Understanding of the tiltrotor configuration has grown significantly but more remains to 
be done. One of the venues where additional research is needed is soft-inplane hingeless tiltrotors. 
An analysis tool is needed to accurately predict the instabilities of ground and air resonance. 
Once this analysis is in place, numerous studies need to be conducted to examine the effects of 
various parameters on these instabilities. One possible source for a solution comes from aeroelastic 
couplings: composite couplings in the wing and kinematic couplings in the blades. There is a clear 
place for this research to add to the body of knowledge. 
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Tiltrotor Modeling 
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3.1 Approach Overview 


A semi-span tiltrotor is modeled using the right wing, pylon, and rotor assembly. The root 
of the wing is assumed to be fixed in space as for a wind tunnel model. The wing is modeled as an 
elastic finite element beam. Three displacements in the x, y, and z directions are modeled between 
the pylon center of gravity and the hub and three additional displacements are allowed between 
the pylon center of gravity and the wing attachment point. The pylon is rigid and attached to the 
wing with a downstop spring. The blades are assumed to be rigid with flap and lag degrees of 
freedom. Any desired twist can be used. The model is pictured in Fig. 1 . Response calculations 
are performed to obtain the steady coning and inflow as well as one of ct and 9 75 the other of 
which is assumed to be known. An eigenvalue analysis provides stability information. 

The model itself is formulated using Hamilton’s principle to derive the equations of the 
motion for the entire model. 


611 = 



-6T- 6W) dt = 0 


( 1 ) 


Since energies are additive, the wing, pylon, and blades can be treated separately. The kinetic, 
strain, and work energies of the blades with arbitrary hub motion are derived initially. The motion 
of elastic blades with hub motion is modeled following the method described in the UMARC 
manual [3] . The hub can be gimballed, hingeless, or articulated. A rigid blade simplification 
is made: the kinetic energy and work energy terms can then be integrated out the length of the 
blade. A simple strain energy formulation is used including springs for the flap and lag degrees of 
freedom as well as the gimbal. The wing is modeled using the finite element method as a simple 
beam with five degrees of freedom at each node — extension in the wing is constrained to zero. 
The rigid pylon serves to connect the hub motion with the wing-tip motion. The contributions to 
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the kinetic, strain, and work energies from the blades are then added to the contributions from the 
pylon and wing. Hamilton's principle is applied to produce global equations of motion. 

3.2 Blade and Hub 

The blade is modeled originally as an elastic beam. Any displacement and rotation is 
allowed, each a function of position on the blade. The degrees of freedom used are given in the 
vector 

[ u e v w <t> x F y F z F <p s a s ip s J (2) 

where xp, yp, and zp are hub displacements, (p s , a s , and ib s are hub rotations, and u e , u, and w 
are displacements in the cross-sectional frame as shown in Fig. 2. Hamilton’s principle is used to 
derive the equations of motion for the blade with six degrees of freedom for the hub motion. 

Note that throughout the derivation whenever allowable, an ordering scheme was applied. 
The orders of the various variables as given in Appendix 8.3 were adapted from the work of Smith 
and Nixon [44] , [35] . The resulting equations of motion for the blade with hub motion were of 
order e\ that is, the largest terms in each of the equations were of order e. It was decided to retain 
all terms of orders e , 52 , e 2 , and e 2 and only linear terms of order e 3 . This convention is applied 
throughout the derivation of the blade and hub matrices. 

3.2.1 Kinetic Energy 

The variation of kinetic energy is given by 

6T = j R Jj pV ■ 6Vdr)d( dx (3) 

To find the velocity V, the position vector is determined and differentiated with respect to time. 
This position vector is expressed as the displacement from a fixed point in space to any point on 
the undeformed elastic axis as expressed in the inertial coordinate system. This vector r is given 
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by the following equation. 

x F 


r = 


vf 
ZF 


1 T 


+ [TwiY { T hw f [T r hY [Tur] 1 [T a 


u 


v cos &o + w sin 9q 
— v sin 00 + iv cos 0 q 



The transformation matrices used in this formula are given explicitly in Appendix 8.4. In general 
[T w i] accounts for the rotation of the hub given by <f> s , a s , and ip s . [Th w ] accounts for the pylon 
rotation a p . \T r h] accounts for the azimuthal rotation ip. [T ur ] accounts for precone and gimbal 
angle as discussed below. [T cu ] accounts for the collective pitch setting and [ Td c ] expresses the 
displacement due to the elastic bending of the blade, rf and £ are cross-sectional coordinates and x 
is the distance along the blade. Foreshortening effects are included so that the degree of freedom 
u is the expressed as the difference between the elastic axial deflection u e and the kinematic axial 
deflection up. 

u = u e — up — u e — I (v f Sv r + w'6w f ) dx ( 5 ) 

Jo 

Because this position vector is expressed in the inertial frame, the velocity V can be calculated 
without recourse to differentiating unit vectors for the frame of reference. 

The variation of velocity is then calculated. Since the velocity includes time derivatives 
of the degrees of freedom, the variation of velocity SV includes terms such as 6v c . Integration by 
parts is used to eliminate such terms. 

Jt ! I j p {T Vc 6v c + Ty c 8v c ) dA dx dt = a: ip P'-Ft 7 p^iAtedt (6) 

A A 

This integration by parts allows the kinetic energy to be expressed as 


£ 8T dt — J m (Tp + ^ T q .6qi j dx dt 


(7) 


where q t ranges over all the degrees of freedom and their derivatives with respect to blade position 


x and Tp are the foreshortening terms. 



3.2.2 Aerodynamics 
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The virtual work SW needed for Hamilton’s principle comes from the aerodynamic forcing 
on the blades assumed here to be quasi-steady. Uniform inflow is also assumed. The aerodynamic 
forces at an arbitrary blade station can be expressed as lift, drag, and pitching moment. 

L = i pV 2 cC t 

D = i P V 2 cC d (8) 

M = \ P V 2 c 2 Cm 

L w , L v , and L u are the external loads acting along the deformed axis. M $ is the moment about 
the deformed elastic axis. These loads are shown in Fig. 3. Assume forms for the aerodynamic 
coefficients as 


Cl 

= Co +aa 


c d 

= C d + d,20t 2 

(9) 

C m 

= C mac + f\(%- 



The first order polynomial term d\ \a\ in the expression for the coefficient for drag is not included. 
An attached flow assumption is made along with a small angle approximation for a so that 


sin a = a: 


cos a = 1 


V = U T 


sin A 
a 


Ur 

U t 

_Up 

U T 


( 10 ) 


These substitutions along with non-dimensionalization give the following forms for the 
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aerodynamic forces: 

1 7 A 

Lu - (-CdUrUfi) ( 11 ) 

L v = ^ihUp 4 - — (— CdUj* — d,2Up — cqUtUp) 

1 1 'y Ik 

= —Q^hUrUp + (coUt — CdUrUp) 

l 1 /y/. 

^0 = ^hedUrUp 4- — (cCmacAj- — e^coAj- — cfiUpUp 4- eaCdUrUp) 

The incident velocity comes from two major sources: the wind velocity and the blade 
velocity. The position on the blade with respect to the hub can be calculated as in the previous 

section. This position is differentiated with respect to time to give the blade velocities. 

u — (y* + w' (77 cos 9 1 — £ sin ) 

— (li)' — v'cfrj (77 sin 0i 4- C cos0i) 
v — cj> (77 sin 0i A £ cos 0i ) 
w 4- <j> (77 cos 0i — Q sin0i) 
where 0 i is the sum of the elastic twist and the pitch setting at this blade position, 9i = e 0 + <t>. 

The velocity of the wind comes from two sources: inflow and forward velocity. The inflow term 

A is vertical in the hub plane. The forward velocity component is horizontal in the inertial frame. 




Note that this equates to a different definition of /x than is used in the helicopter community: 

A 4 — f° r helicopters where a s is the tilt of the shaft, but /x = ^ for tiltrotors. V w as shown 



“strip” analysis based on the angle of attack at the section three-quarter chord location; rj T is the 



24 

chordwise offset between the elastic axis and the 3/4 chord position. Up , Up, and Up are found by 
totaling the velocities above and transforming them to the deformed frame. 

I Ut } = [T dc \ [T cu ] |^|= [T dc \ [T cu ] (V bh - V w ) (15) 

Once these velocities are expressed in the deformed frame, they can be substituted into the forms 
for the aerodynamic forces as given above to give the virtual work terms. 

3.2.3 Gimbal Modeling 

To give some degree of confidence that this tiltrotor model is realistic, it was necessary 
to model a gimballed hub. The vast majority of wind tunnel data and other analyses for tiltrotors 
are based on gimballed hubs since the tiltrotors flying today are use gimballed rotors. The gimbal 
is modeled throughout following the method described in [35] . The gimbal is treated as two 
additional degrees of freedom of hub motion where 

@G = 0GC cos T- PgS s * n V* (16) 

This (3 C is added into the transformation matrices with precone so that the transformation matrix 

between rotating and undeformed coordinate systems takes the following form. 

cos j3 p + sin (3 p {(3 GC cos ip + j3 GS sin ip) 0 — cos (3 {(5 GG cos ip 4- (5 GS sin ip) + sin /?_. 
[Tur]= 0 10 p 

. cos ftp (Pec cos V’ + Pgs sin xb) — sin d p 0 cos (3 p + sin (3 p (ft GG cos ip + p G g sin ip) 

(17) 

The effect of the gimbal motion is shown in Fig. 4. 

3.2.4 Rigid, Twisted Blade 

The resulting kinetic energy terms and work energy terms are simplified to rigid blade 
coordinates. Many studies have shown that rigid blades can provide a good approximation for 
stability analyses. This simplification allows comparison with other rigid blade analyses and with 
equations of motion derived by hand. Also, since the present version of this analysis does not 
include modal analysis, using a rigid blade keeps the number of degrees of freedom to a minimum. 
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The primary elastic blade degrees of freedom are simplified to rigid flap and lag degrees 
of freedom. Extension and elastic twist are not modeled in the rigid blade code. 


U e 

= 4> = o 


Vb 

oT 

1 

.2 

*53 

1 

II 

(18) 

W b 

= sin (3(x — ef) 



All of the elastic blade terms are simplified in a similar manner; the substitutions are explicitly 
given in Appendix 8.5. The resulting degrees of freedom for the blade / hub system in the rotating 
frame are 

[ [3 (, x F xjf Z F a s 4> s ip s (3 gc 0 gs J (19) 

The chordwise offset between the blade center of gravity and the elastic axis and the flapwise and 
chordwise mass moments of inertia per unit length along the blade are preserved. The effects of the 
foreshortening of the blade are also retained. The rigid blade assumption also allows the integrals 
associated with this foreshortening to be evaluated. For example, 

j j v v dr] dQ — x sin 2 Q (20) 

'A 

The results of the other integrals associated with foreshortening are given in Appendix 8.5. 

Since the blade is assumed to be rigid, it is then possible to integrate the work and kinetic 
energy terms out the length of the blade. In both cases the terms are explicit in x\ the implicit 
variation with x comes with blade twist. Both the kinetic energy terms and the work terms are 
functions of the pitch at station x out the length of the blade. Since the blade can be twisted, 

$o = #o (^)- One choice would be to assume a polynomial expansion or similar known form for 
do ( x ). However, in this case no assumption is made about #o- The blade twist do is set equal to d c 
at the blade root plus an arbitrary function of x d tw ( x ) . Since do is not assumed to be a small angle, 
do appears throughout as cos do — cos ( d c 4- d tw ) and sin do = sin ( d c + d tw ). For an untwisted 
blade approximation it is trivial to set 6t w to zero throughout. Each individual term is expressed 



as c / mx nl cos" 2 0 tw sin" 3 9 tw dx where c is not a function of x and the constants nl, n2, and n3 
vary between 0 and 2. These integrals are left unevaluated until a twist function is specified with 
the other rotor parameters. This method actually results in fewer terms than assuming a linear twist 
and integrating into a closed form solution; the need to integrate numerically out the blade is also 
removed. 

3.2.5 Strain Energy 

Since a rigid blade assumption has been made, the strain energy for the blade is limited to 
the springs at the root of the blade. 

SU = k/3 (/? — /3 p ) 8(3 + A :q £ 8(( 

+ k P a Pgc 60 gc cos2 $ 

+ { k 0 (: 0GC 6 Pgs + k p G Pgs 6 Pgc) cost/) sin i/> 

+ k 0 a Pas bPas sin 2 ip (21) 

3.2.6 Structural Coupling 

Structural coupling between the blade flap and lag is modeled as shown in Fig. 5. This 
allows a better approximation of a hingeless blade as R is allowed to vary between 0 and 1. The 
structural coupling terms appear only in the stiffness matrix in the flap and lag equations [2] . 

3.2.7 Linearization 

These strain, kinetic, and work energies are linearized by assuming that each of the 
degrees of freedom is equal to a steady value plus a perturbation: for example, /? = /? + 60. This 
substitution is exact-no simplification is assumed here. The next step, however, is to assume that 
products of perturbations are small: 50 ■ = 0, for example. As long as the response routine 

is of sufficient accuracy, this assumption is reasonable. The resulting steady equations become a 
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response routine. If the hub were assumed fixed in space, Hamilton’s principle could be applied 
at this point to produce equations of motion. With pylon and wing motion, the kinetic, strain, and 
work energies from the blade and hub must be added into the kinetic, strain, and work energies 
for the pylon and wing before Hamilton’s principle can be applied. However, it is convenient to 
express the linearized perturbation values for the blade/hub kinetic, strain, and work energies in the 
form of preliminary [ Mbiade ]. [C blade], and [K blade] matrices. A Fourier coordinate transformation 
is applied to transform the rotating system coordinates for each blade into fixed-frame coordinates. 
The resulting fixed-frame matrices are assembled with the resulting wing and pylon kinetic, strain, 
and work energies to give equations of motion for the entire model as described below. 

[M] {<?} + [<7] {<?} + [I<\ { 9 } = {F} (22) 

3.3 Pylon and Wing Attachment 

The pylon is modeled as a rigid connection between the hub and the wing. The kinetic 
energy associated with the motion of the pylon mass and inertia is added to the kinetic energy 
associated with the wing-tip motion and the hub motion resulting in additional terms in the total 
mass matrix. 

ip ] [ m Fx 0 0 0 0 0 ' 

VP 0 m Fy 0 0 0 0 

ip l 0 0 m Fz 0 0 0 

'4> P [ 0 0 0 /^00 

a P 0 0 0 0 ’ 7 a , 0 

' i} P J [0 0 0 00 " 

where the degrees of freedom xptoipp are related to the wing-tip degrees of freedom as described 
below. To accurately predict the coupling between the wing torsion and vertical bending, three 
displacements in the x, y, and 2 directions are allowed between the pylon center of gravity and the 
hub (xcg, Vcg, h) and three additional displacements between the pylon center of gravity and the 
wing attachment point (x pcg , ypcg, z pcg )(see Fig. I above). The offsets in the 2 direction h and z ix:g 
are important because the rotor and pylon mass offset from the wing couples wing bending and 
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torsion. These offsets are taken into account in the assembly process. 

The pylon is attached to the wing with a downstop spring, a torsional spring between the 
wing-tip and the pylon in the direction of the pylon rotation. As the pylon is rotated forward into 
forward flight mode on current tiltrotors, the pylon is locked into position with the wing. This 
spring and locking mechanism are modeled by allowing motion between the pylon pitch and the 
wing-tip pitch thus creating an additional degree of freedom for the total system. The locking 
mechanism is modeled by allowing the spring to be very large in forward flight. The hub motion 
is calculated in terms of the pylon pitch <f>N+2 while the wing motion is calculated in terms of the 
wing-tip pitch <£^ +1 . These two degrees of freedom are related by the strain energy associated 
with the downstop spring. This strain energy adds into the strain energy for the blades and the wing 
providing additional terms in the total stiffness matrix. 



3.4 Elastic Wing 

The wing is assumed to be an elastic beam and modelled using the finite element method 
as shown in Fig. 6. Each node has 5 degrees of freedom: extension in the wing is constrained to 
zero. At each node the degrees of freedom are 

{<7 node} = [ V V W W <j> j (25) 

in the directions shown above. 

3.4.1 Kinetic Energy 

The kinetic energy for the beam element is, by definition, 

ST et = / pSu T iidV (26) 

Jv 

For a beam with an offset ycM between the center of gravity and the elastic axis, the displacement 
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vector is given as 


u = 


u 

V 


w + vcm4> 

The shape functions used to model the wing are as follows: 
n(x,t) = [ (1 -f) ?]{“*} 

= 1 f 1 -*) f]{^ } 


v{x,t) = { 


w(x,t) = 


1 - 3(f) 2 4-2 (t) 3 

i((f)-2(ff+(f) 3 ) 
Hi) 2 -Hi ) 3 
t (- (i ) 2 + (i) 3 ) 

i-Hif+Hi ) 3 
(ti) - Hi ) 2 + (if) 

, 3 (x) 2 - 2 (r) 3 

l i(-(f) 2 + (i) 3 ) 




Z? 



( 27 ) 


( 28 ) 


Making these substitutions and carrying out the math gives an element mass matrix. The complete 
wing matrices are given in Appendix 8.6. Note that the wing vertical bending and torsion are 


coupled due to the chordwise offset between the wing center of gravity and the elastic axis. It 
is assumed that any offset between the center of gravity and the elastic axis in the 2 direction 


(perpendicular to the chord) is small enough to be neglected. 




/ 

u A 

6v a 


v A 

6v a 
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8w a 
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8ub 

> [M 5 ] < 

&A > 
Ub 

8v b 
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8v' b 
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8w b 


W B 

6w' b 


U)' B 

< 8<Pb * 


, <t>B J 


= SU T [M 5 ] U 


( 29 ) 


Though the extension degree of freedom u is shown in the degree of freedom vector here, the 



extension in the wing is later constrained to be zero. 
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3.4.2 Strain Energy 

The definition of strain energy can be applied with the shape functions given above to 
calculate an element stiffness matrix. 



6U el = U T [K s ]U (31) 


This form of the strain energy for the wing is modified with the inclusion of the aeroelastic 
couplings as described below. 

3.4.3 Work Energy 


The angle of attack of the wing is given by 


w 


a = t*o + 4> cos A — — — w' sin A 


(32) 


Note that wing sweep A is included which is important because Nixon found wing sweep to be 
destabilizing [34] . 

The forces are given by 


P= < 


( I 1 

*->u 


' 0 

L v 


0 

Mz 

► = < 

0 

1 L w 

jpV 2 cci n (ao + ^cos A — — u/sin A) 

My 


0 

M* J 


k \pV 2 cci a e (ao + 0cos A — y — ti/sin A) 


(33) 


where e is the chordwise distance of the elastic axis behind the aerodynamic center 
The virtual work for the element is given by definition as 

8W el = [ PSUdx 

Jo 

Substituting in the same shape functions as above gives aerodynamic stiffness and damping 


( 34 ) 
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matrices for the wing element. These matrices are given explicitly in Appendix 8.6. 

6W el = 8U ({F 4 } + [I <- 4 ] U + [C A ] U^j (35) 

3.5 Assembly 

Once the kinetic, strain, and work energies are calculated from the blade, hub, pylon, 
and wing, Hamilton’s principle can be applied to assemble all of these energies into total system 
equations of motion. The [Mbiade], [Cbiade\, arid [K blade) from the rotor are combined into total 
system [A/], [C], and [K] by equating the wing-tip degrees of freedom to the hub motion degrees 
of freedom. Since the pylon is rigid, this is just a one-to-one correspondence that can be expressed 
as a transformation matrix. The single exception is due to the addition of the downstop spring 
which is modeled as an extra torsion degree of freedom in the total degree of freedom vector. 

Since the hub degrees of freedom xp, yp, zp , a 3 , (j> s , and ip s from the hub and blade 
equations are defined in the inertial frame, the hub degrees of freedom can be equated to the wing 
tip degrees of freedom in a one-to-one fashion. 


Xh 


' 0 

1 

0 

0 

0 

0 ' 


/ > 

u 


/ 

u 

Vh 


1 

0 

0 

0 

0 

0 


V 


V 

Z h 


0 

0 

0 

1 

0 

0 


v' 


v f 

—<*h 

> — 

0 

0 

0 

0 

0 

1 

< 

IV 

> = [EE] < 

w 

~<Ph 


0 

0 

0 

0 

1 

0 


w’ 


iv r 

1>h j 


0 

0 

-1 

0 

0 

0 


<$> 


0 J 


This transformation becomes complicated by the six offsets. A rotation of the wing tip produces 
displacements at the hub because of the offsets xqg, VCG , and h between the pylon center of 
gravity and the hub and because of the offsets x pcg , y pcg , and z pcg between the pylon center of 
gravity and the wing tip. The exact transformations are given in Appendix 8.7. 

A numerical Fourier Coordinate Transformation is applied to the assembled mass, stiffness, 
and damping matrices to transform the rotating frame flap and lag into fixed frame coordinates 
[8] . The resulting total degree of freedom vector consists of the five degrees of freedom at each 
wing finite element node, the blade flap and lag terms, and the gimbal terms. Note the presence of 
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0/v+2 * n the last blade element which is the extra 4 > degree of freedom due to the presence of the 
downstop spring. 

{ qtot } = [ Eh ... EIn blade El ^ + 1 gimbal J T (37) 


[Eh \ = 

[ V\ v\ W\ iu\ (f> J 


[ EIn J = 

[ 0N J 


[ EIn+ 1 J = 

[ V N+ 1 v N+l w N+\ w N+l ^N+l 

&N+2 

[blade J = 

[00 01c 013 Co Cl c Cls •• 

• J 


[ gimbal J = [ 0 GC l3 GS J 

3.6 Response 

Once the total system matrices are assembled, the response is needed, in this case the 
coning and pitch setting. The steady values of the blade and hub degrees of freedom are assumed to 
be zero in all cases except coning. A value for #75 or c? must be given. The inflow equation comes 
from momentum theory while the coning and thrust equations are found from the steady terms in 
the equations of motion as described above. The coupled response routine uses a Newton-Raphson 
iteration to calculate the coning ( 3 0 , the resulting inflow A, and the unknown of #75 or cr- 

3.7 Stability Analysis 

The mass, damping, and stiffness matrices for the entire system can now be calculated. The 
forcing vector is set to zero for the stability calculation and the resulting matrices are transformed 
to state space form. 

m {ii> + [C] {*> + [K] M = {0} =1 [ M j { 9 } + [ - M j {„} = {0} (38) 
The stability of the degrees of freedom is calculated from the eigenvalues. Currently stability 
analyses can be conducted by varying either the rotor rotational speed. Cl, or the forward speed of 



the tiltrotor, V. 
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3,8 Modeling of Couplings 


The sections above describe the modeling and analysis approach for a baseline model. 
This section describes how the couplings in the wing and blade are implemented. 

The composite couplings in the wing modify the wing strain energy as shown below 


U e i — g j EA ^ + GJ {^p j 4- EI Z ^ -f EI y ^ w ^ + 2 (^K^v + K^w j <p 


dx 
( 39 ) 


Thus A45 couples wing vertical bending and twist while K40 couples wing chordwise bending and 
twist. The sign convention used here is that positive wing vertical bending-twist coupling is when 
the blade bends up and the nose of the airfoil pitches up (^45 > 0 ). Wing chord bending-twist 
coupling is defined as positive when the blade bends forward and the nose pitches up (/Qg > 0). 

The kinematic couplings in the blades are included as changes in the effective twist. 

#0 = 0 q ~ p tan 63 — ( tan £4 + tan 8 z G (P GG cos xp -b p GS sin ip) 


= 0 o — pk P0 — <^k Pc -b k p p G (p GG cos xp -b / 3 G g sin xp) ( 40 ) 

These coupling terms appear everywhere that the pitch does. For example, the transformation 

matrix between the cross-section frame of reference and the undeformed frame is given by 

‘10 0 

[Tcu] = 0 cos 0 O sin ( 41 ) 

0 — sin 0o cos 0 q 

With the blade couplings and gimbal degrees of freedom included, the [T cu ] matrix becomes 
10 0 

T r _u=\ 0 (1 -l-T ( i u )cos6>o + r c 2 u sin6'o -7^cos0 o + (l + T^) sin<9 0 | (42) 

0 cos 6q- (1 + 7^) sin 6 0 ( 1 + T^ u ) cos 6 0 -f sin 8 0 


where 


Tea — hpfi (3 kpQ £ + (kpp (3 4 - kpQ £) (k p /3 G ( 3 gg cos xp + kpp G p G $ sin xp) 
Tcu ~ kp /3 (3 -b kp^ £ — kpp G (3 gg cos xp — kpp G p G $ sin xp 
+k pf 3 (3 kp r £ {kppe Pec cost/; -b k p/3G p GS sint p) 


( 43 ) 
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The transformation matrix between the deformed frame and the cross-section frame is also altered 


by the inclusion of the blade couplings. This matrix is given explicitly in Appendix 8 . 4 . The 63 
terms couple blade flap and pitch while 64 couples blade lag and pitch. Positive pitch-flap coupling 
is defined to be when the blade flaps up and the nose pitches down (£3 > 0 ). Positive pitch-lag 
coupling is when the blade lags back and the nose pitches down (64 < 0 ). The blade couplings are 
included in the response calculations as well as in the stability calculations. 



4 Validation 
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This section details the validation study which was conducted to give confidence in the 
present analysis. Comparisons were made for both a Boeing model which was hingeless and 
soft-inplane and for a Bell model which was gimballed and stiff-inplane. Several abbreviations 
are used in labeling the modes: WV=wing vertical bending, WC=wing chord bending, WT=wing 
torsion, LFL=low frequency rotor lag, HFL=high frequency rotor lag, CL=collective lag, LFF=low 
frequency rotor flap, HFF=high frequency rotor flap, CF=collective rotor flap, LFG=low 
frequency gimbal flapping motion, and HFG=high frequency gimbal flapping motion. 

4.1 Comparison with Hingeless Rotor for Air Resonance 

In 1973 Boeing conducted two wind-tunnel tests of the Model 222 rotor in the NASA 
Ames 40’ by 80' wind tunnel [26] . This was a full-scale, semi-span model with a 25 -foot-diameter 
rotor (Fig. 7). The main purpose of this tunnel entry was to gather performance and aeroelastic data 
about the Model 222 to validate Boeing's analysis of that rotor. Regions of degraded stability in 
forward flight were observed. Testing was done only in forward flight mode; no ground resonance 
data was accumulated. 

Four different wind tunnel speeds were used: 50 knots, 100 knots, 140 knots, and 192 
knots. Damping for the wing modes and frequencies for the wing vertical bending were gathered. 
The data from these wind tunnel tests along with the Boeing analysis of that time are the primary 
sources of validation for the current analysis. More data was accumulated for the 100 knot run so 
that will be the focus of the validation and preliminary studies. 

4.1.1 Fifty Knots 

Fig. 8 shows the damping for the wing modes for the Boeing damping as measured at 50 
knots. The frequencies are not reported for 50 knots. The circles indicate data gathered at 50 knots 
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while the triangles show data gathered at 60 knots in the wind tunnel. Both analyses shown are for 
50 knots wind tunnel speed. The wing vertical bending damping shown by the current analysis 
matches Boeing’s analysis and the wind tunnel data well. The wing chord damping analyses differ 
by approximately 2% critical damping offset. The wing torsion mode is also predicted to be less 
damped by the current analysis. The high frequency lag mode is also shown here for the current 
analysis because of the high degree of coupling between this mode and the torsion mode around 
280 RPM. 

4.1.2 One-Hundred Knots 

Shown in Fig. 9 and 10 are the frequencies and damping as given for the wind tunnel test, 
the Boeing analysis, and the current analysis at 100 knots [26] . The circles are the wind tunnel data 
that was generated. Boeing’s analysis is shown in dotted lines while the current analysis is solid. 
The frequencies predicted by Boeing’s analysis and the current analysis match well. In general the 
wing frequencies predicted by the current analysis are higher than the wing frequencies predicted 
by Boeing’s analysis. One possible source of error here comes from the lack of data given in the 
report [26] : the wing frequencies (blades off) are reported as 2.5 Hz for wing vertical bending, 4.5 
Hz for wing chord bending, and 11.3 Hz for wing torsion. Milues for the stiffnesses EI y , EI Z , and 
GJ are not explicitly given. These parameters were estimated based on the weight of the wing and 
pylon. It is assumed that Boeing’s analysis is based on the actual stiffnesses of the wing design. 

One discrepancy is in the wing torsion mode. Boeing’s wing torsion mode is shown to 
vary from 4 Hz at 100 RPM to more than 10 Hz at 600 RPM. A body frequency varying by that 
magnitude with change in a rotor rotational speed is considered unlikely. However, this mode is 
matched well by the current analysis’ high frequency lag mode for lower RPM and is matched less 
well at high RPM by the current analysis’ wing torsion mode. The vertical lines at 260 RPM and 
380 RPM shown in the current analysis are regions where the wing torsion mode couples very 
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strongly with high frequency flap around 260 RPM and with high frequency lag around 380 RPM. 
The high frequency flap and lag mode frequencies clearly increase with rotor RPM. In the regions 
indicated, choosing one frequency for high frequency flap and the other frequency for wing torsion 
is inaccurate. Barkai and Rand label both frequencies in these regions with both modes, e.g., high 
frequency flap/wing torsion mode, which is one way of indicating these strongly coupled modes 
[ 1 ]. 

The damping results show more discrepancies with the Boeing analysis than do the 
frequency results. As in the 50 knot case, the high frequency lag mode is also shown for the current 
analysis. Since this stability data is shown as a critical damping value, the discrepancy discussed 
for the wing torsion mode frequency above will affect the wing torsion damping shown here. The 
wing vertical damping predicted matches the wind tunnel data fairly well and certainly shows a 
region of air resonance as seen in the wind tunnel. The shift of the instability to higher RPM is 
probably due to the higher frequency predicted for the wing vertical bending mode. The wing 
chord mode damping as predicted by the current analysis is similar in tendency to Boeing’s analysis 
but differs in magnitude by approximately 2%. The wing torsion mode is similarly affected. This 
shift in damping is currently not understood. 

4.1.3 One-Hundred Forty Knots 

Fig. 11 and 12 show the frequency and damping for the wind tunnel test at 140 knots along 
with the Boeing analysis and the current analysis. The circles indicate wind tunnel data gathered at 
140 knots. This shows that the wing vertical bending damping correlates well with the wind tunnel 
data for various forward speeds. Again the frequency data is excellent. The wing chord damping 
analyses differ by the same 2% critical damping offset at high RPM. The wing torsion mode is also 
predicted to be less damped by the current analysis. The high frequency lag mode is also shown 
here for the current analysis because of the high degree of coupling between this mode and the 



torsion mode around 280 RPM. 
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4.1.4 Discussion 

The calculation of the steady coning, inflow, and pitch or thrust requires that one of the 
pitch and thrust be given over the rotor rotational speeds to be tested. The Boeing report included 
the pitch settings for the wind tunnel testing but not the pitch settings used in the accompanying 
analysis. Data for the 6 75 setting was only available for 7 test points between 320 and 412 knots 
for the 50 knot test, for 24 test points between 104 and 470 knots for the 100 knot test, for 11 test 
points between 384 and 491 knots for the 140 knot case, and for 7 test points between 315 and 386 
knots for the 192 knot test. Linear interpolation was used to estimate #75 values outside of these 
ranges. Because of the additional data available for the 100 knot case, all of the results shown in 
forward flight mode in the next section were run at 100 knots. This data is shown graphically in 
Fig. 13 where the crosses mark values given and the line shows data used for this analysis. 

Choices of which modes to model also affect the results. For example, Fig 14. shows 
how including the collective lag mode in the stability analysis changes the behavior of the wing 
chord damping. The dotted lines are the wing chord bending and rotor collective lag modes as 
predicted by the current analysis for frequencies and damping. The solid lines are taken from the 
current analysis results having excluded the collective lag mode. This change in behavior calls 
into question exactly what modeling assumptions were made for the Boeing analysis, specifically 
whether or not the high frequency rotor flap and lag were included. This change in behavior of the 
wing chord mode is not expected from helicopter theory: the collective lag mode does not couple 
with the body modes to create an air resonance situation as is seen here. Only the low frequency 
lag mode is crucial for determining the stability of a helicopter from air and ground resonance. 

The airfoil data is not available. Nor are the blade effective hinge offsets or blade 
non-rotating natural frequencies. These parameters have been approximated for use in the current 
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analysis. The damping is quite sensitive to some of the approximated parameters. Lack of data 
such as this could account for some of the discrepancies between the Boeing analysis and the 
current analysis. 

4.2 Comparison with Gimballed Rotor for Whirl Flutter Analysis 

Confidence with a model’s accuracy comes when the model predicts the same behavior 
as seen in wind tunnel results and other analyses. Since the majority of research conducted 
with tiltrotors concerns whirl flutter, it is important that the current analysis predict whirl flutter 
accurately even though this is not the primary emphasis of this research. The current model 
includes the gimballed hub specifically for this reason. The gimballed hub modeled here is a Bell 
Helicopter hub on the same NASA wing as used above. The first comparison is made with the 
response predictions from Nixon’s elastic blade code in UMARC as in Fig. 15 [35] , [3] . The c ? 
values are given and the pitch settings are calculated. This result gives confidence that the response 
calculation is consistent with that of other analyses. 

The results to be compared next are from Nixon’s rigid blade model [35] . Fig. 16 shows 
the frequencies for the whirl flutter analysis which is plotted with increasing velocity rather than 
increasing rotor rotational speed as for ground resonance. The solid lines are the current analysis 
while the dotted lines represent Nixon’s analysis. The collective lag motion is not shown. The 
agreement between the prediction of the frequencies is excellent. The corresponding damping 
values are given in Fig. 17. Again, the current analysis is shown with solid lines while Nixon’s 
analysis is shown using dotted lines. The agreement in the damping is not good. The wing mode 
damping is smaller for the current analysis prediction as was true for the soft-inplane validation 
study. In addition, this analysis predicts more lag damping and less flap damping than does Nixon’s 


analysis. 



4,3 Summary of Validation 
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The air resonance trends shown by the Boeing wind tunnel test from 1973 are well 
predicted. The wind tunnel data that was obtained is matched well. The differences in modeling 
coupled with the lack of data could easily account for the differences seen in the Boeing analysis 
and the current analysis. The response routine matches very well with results from Johnson and 
Nixon [17] , [35] . The whirl flutter frequencies are also in very good agreement with Nixon's 
analysis. The damping values for whirl flutter remain suspect and require additional attention. 
These validation studies, though not totally completed, give enough confidence in the prediction of 
tiltrotor air resonance to justify a first parametric study to examine the use of couplings to stabilize 
the air resonance instabilities shown. 



5 Preliminary Results and Discussion 
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The preliminary results presented here can be divided into four parts based on the couplings 
discussed: wing vertical bending-twist coupling, wing chordwise bending-twist coupling, blade 
pitch-flap coupling, and blade pitch-lag coupling. For each coupling the situations of air resonance 
in forward flight mode, air resonance in hover mode, and ground resonance in a non-thrusting 
case must be considered. Unlike in a helicopter where the zero-thrust case also entails adding 
the landing gear stiffness and damping if present, the ground contact case for the present model 
tiltrotor does not change the stiffness or damping matrices. The hover case is for a nominal value 
of = 0.2 for all RPMs. This gives values of 6 75 of around 15 degrees. 

The baseline frequency and damping for the 100 knots case is given above in the validation 
section. The baseline frequencies and damping for the hover case are given in Fig. 19 and 20. In 
the forward flight mode the wing vertical bending is very lightly damped. However, in the hover 
mode, the wing vertical bending couples with the vertical plunge motion of the rotor which is very 
heavily damped. This increase in damping in the wing vertical mode eliminates the air resonance 
seen in forward flight mode. Note also that the nominal RPM for the hover configuration changes 
to 551 RPM from the forward flight nominal RPM of 386. This hover case assumes that the wind 
tunnel velocity is zero since the analysis currently assumes constant-coefficient matrices. 

The baseline frequencies and damping for the zero-thrust case are given in Fig. 21 and 22. 
The zero-thrust case is also in helicopter mode so again the nominal RPM is 551 RPM. Note that 
the cyclic lag mode is included in these figures. 

For each type of coupling, the first question examined was whether the positive and 
negative coupling of that type stabilized, destabilized, or had no effect on air resonance. Couplings 
which destabilized air resonance were no longer considered. If a positive or negative coupling 
was shown to stabilize air resonance (or to have only little effect), the effect of that coupling 
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was examined on hover and zero-thrust conditions. If hopeful conditions were found for a given 
coupling for forward flight, hover, and zero thrust, then whirl flutter was considered. 

5.1 Wing Vertical Bending-Twist Coupling 

Wing vertical bending-twist coupling was defined as positive with bending up coupled to 
nose-up twist, i.e., > 0. Negative beamwise bending-twist coupling < 0) destabilizes air 

resonance. Positive beamwise bending-twist coupling (#45 > 0) tends to stabilize air resonance 
(Fig. 23). This destabilizes the zero-thrust and hover conditions slightly but since this tiltrotor is 
extremely stable in these conditions the destabilizing trend is unimportant (Fig. 24 and 25). 

5.2 Wing Chordwise Bending-Twist Coupling 

Recall that the wing chordwise bending-twist coupling was defined as positive when the 
wing bends forward and twists nose up, or > 0. The effects of chordwise bending-twist 
coupling on air resonance, hover, and zero-thrust conditions was quite small. This coupling did not 
stabilize the existing air resonance in the Boeing baseline model. If this coupling was required for 
some other purpose, it is good to know that this coupling does not significantly harm air or ground 
resonance but neither does it stabilize these conditions. 

5.3 Blade Pitch-Flap Coupling 

Blade pitch-flap coupling is defined above to be positive when blade flapping up 
couples with blade twist nose-down, 63 > 0. Both the positive and negative pitch-flap couplings 
destabilized the air resonance in the baseline case, though the effect was not very strong. 

5.4 Blade Pitch-Lag Coupling 

Let pitch-lag be positive when the blade lags back and pitches nose-down, £4 < 0. 

The negative pitch-lag coupling destabilized the air resonance. However, the positive pitch-lag 



coupling stabilized the air resonance in the baseline rotor (Fig. 26). These graphs show the blade 
low frequency lag and the wing vertical bending. The solid line on these graphs is the baseline 
case, the dot-dashed line indicates —5 degrees and the dotted line shows —10 degrees. The 
pitch-lag coupling on this tiltrotor s air resonance instability is found to be effective. This result is 
in line with what helicopter studies have found [36] . However, this positive pitch-lag coupling in 
helicopters destabilizes the ground resonance instability. 

Fig. 27 and 28 show the hover and zero-thrust damping values for the Boeing model. The 
addition of positive pitch-lag coupling does not destabilize the hover or ground condition. This 
is unlike a the results for a helicopter. One possible cause for this comes from the difference in 
stiffness and damping due to the landing gear on a helicopter which become ineffective in forward 
flight. The wing vertical bending mode modeled here does not show that loss in stiffness or 
damping on take-off. 

5.5 Discussion 

The aeroelastic couplings discussed here impact the stability of the tiltrotor in air resonance. 
The two couplings which show promise for stabilizing air resonance without destabilizing hover 
or ground resonance are positive wing beamwise bending-twist coupling ( A 45 > 0 where the wing 
bends up and pitches nose-up) and positive blade pitch-lag coupling (64 < 0 where the blade 
lags back and pitches nose -down). Before incorporating these couplings into tiltrotors, it must be 
shown that these couplings do not destabilize the tiltrotor in transition or forward flight in hover 
mode nor significantly decrease whirl flutter stability margins. 

Published whirl flutter coupling results show that chordwise bending forward-pitch 
nose down coupling within the wing can add significant stability to the system for whirl flutter 
(^46 < 0). Barkai and Rand showed that wing vertical bending up-pitch nose-up coupling also 



44 

adds to the system stability (/T 45 > 0) [1] . Popelka et al. showed that tailoring the wing so that as 
the wing bends up the pylon pitches nose down could increase the whirl flutter stability (A 45 < 0 ) 
[39] . These effects seem to be altered by the inclusion of blade elasticity which was included for 
Barkai and Rand but not for Popelka et al. 

Nixon showed that flatwise bending-twist coupling in the blades could also increase 
the whirl flutter stability margin when flatwise-bending up coupled with nose-down twist; this 
coupling tends to offset the precone induced pitch-lag coupling which is destabilizing for whirl 
flutter [35] . Tiltrotor blades are highly twisted and use large collective pitches. The chord line 
angle with the rotor plane of rotation is very large. The flatwise bending is defined as bending 
perpendicular to the plane of rotation. Thus, Nixon's study on flatwise-bending twist coupled 
blades indicates that positive pitch-lag coupled blades could detrimentally affect the stability in 
whirl flutter. More research is clearly needed here to assess the effects of these couplings on whirl 


flutter. 



6 Further Research 
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The preliminary results described above indicate a potential for these couplings which 
merits further study. Several refinements are needed and extensive parametric studies should be 
conducted to expand the understanding of if and how these couplings can stabilize air and ground 
resonance for a soft-inplane tiltrotor. The first step is to understand exactly why the gimballed 
rotor damping does not match with Nixon’s predictions. As soon as that is completed, refinements 
can be made to the model itself. 

The first additional feature to be added is the free flight degrees of freedom for the 
fuselage. The model above assumes from the outset that the root of the wing is fixed in space: 
the rigid body motion of the entire tiltrotor is not modeled. These free-flight degrees of freedom 
could also experience air or ground resonance. The study conducted by Srinivas in 1995 included 
these body modes as well as the second rotor and wing. The antisymmetric modes of vibration 
for the wing were therefore modeled as well as the symmetric modes which are modeled here. 
Srinivas showed that for the baseline case he tested, the whirl flutter stability was governed by 
the symmetric modes [46] . It is not clear how the antisymmetric modes will affect ground and air 
resonance of a tiltrotor. These modes need to be addressed and studied. This modification will also 
include horizontal tail effects and free-flight vehicle trim. 

The flexibility in the pylon should also be included. The rigid pylon assumption is not 
unreasonable, but since flexibility in the pylon could change the dynamics of the coupling between 
the rotor and the wing, it is important to model whatever flexibility exists. This is especially 
important for whirl flutter results where the flexibility in the pylon can certainly affect the results. 

One of the serious limitations in the model as it stands now is the need to be either in hover 
or in airplane mode since the mass, damping, and stiffness matrices are assumed to have constant 
coefficients. This limitation eliminates forward flight in helicopter mode and transition. Floquet 
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theory will be implemented to allow non-constant matrices for the stability analysis. The finite 
element method in time will be implemented to calculate the response. Note that the arbitrary pylon 
tilt is already modeled. Other possible improvements would be the use of unsteady aerodynamics. 

The earliest models of ground resonance show that damping must be added in both the 
fixed frame and in the rotating frame to stabilize air and ground resonance. In the case of a 
helicopter this usually entails lag dampers on the blade and dampers on the landing gear. In the 
case of a tiltrotor, damping in the rotating frame could still be provided by lag dampers on the 
blades, but damping in the fixed frame for a tiltrotor could pose more problems since adding 
damping to the lightly damped wing modes is less straightforward. The couplings discussed above 
could certainly help the stability. However if these couplings prove to be inadequate to stabilize 
a soft-inplane tiltrotor, another suggestions is to use other damping in the wing to increase the 
fixed-frame damping. Parametric studies need to be conducted to establish how much damping is 
needed in the fixed frame and how much is already accounted for from aerodynamics. 



7 Proposed Time Line, Publications, and Presentations 
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Validation studies will be completed by the fall semester. The difference in the damping 
values for the wing chord and wing torsion modes must be accounted for before an abstract can 
be submitted. The first abstract deadline of interest would be for the AIAA SDM conference on 
August 1, 1999. The abstract deadline for the AHS Forum 2000 is September 1, 1999. 

With a completely validated model, the fuselage rigid body motion will be added. The 
preliminary studies of the blade couplings 6 s and 64 and the wing couplings A 45 and /v 46 must be 
completed. These studies will have accomplished the tasks of exploring the possibility of air and 
ground resonance in a tiltrotor and evaluating the potential for composite couplings to stabilize 
these instabilities. If the validation and parametric studies are completed, an abstract will be 
submitted for the AHS Forum 2000. 

Once these studies are completed, the additional model refinements will need to be 
implemented. Flexibility in the pylon, Floquet theory in stability and response codes, unsteady 
aerodynamics, and dynamic inflow must all be implemented. The goal is to have these in place 
by summer, 2000. At this point extensive parametric studies must be completed to completely 
understand the phenomena of air and ground resonance in a tiltrotor and to determine whether the 
couplings can stabilize a soft-inplane hingeless tiltrotor. Work will completed on these tasks by the 
end of August, 2000. The dissertation will written and defended in Fall, 2000. 



8 Appendices 


8.1 Figures 



Fig. 1. Trancendental Model 1-G 




Fig. 3. Teetering Hub 



Fig. 4. XV- 1 5 at NASA Dryden Research Center 
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Fig. 8. Frequencies for ground resonance 
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Fig. 9. Damping for ground resonance 
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Fig. 10. Frequency and damping for whirl flutter stability 
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Fig. 11. Tailored wing or blade with box-beam spar 
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Fig. 19. Model Description 



Fig. 21. Aerodynamic Forces in the Deformed Frame 






Fig, 23. Blade Structural Coupling Model 
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Fig. 26. Wing Damping, Boeing Model, 50 knots 



Fig. 27. Frequencies for Boeing Model, 100 knots 
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Fig. 30. Damping for Boeing Model, 100 knots 
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Fig. 34. Frequencies for whirl flutter from Nixon and current analysis 
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Fig. 35. Damping results for whirl flutter from Nixon and current analysis 
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Fig. 36. 9 75 for RPM sweep of Boeing Model 222 in hover, no couplings 
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Fig. 40. Zero-thrust damping ratios, Boeing model, no couplings 
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Fig. 46 . Positive 6 4 coupling for zero thrust, baseline, 64 = —5 degrees, and 64 = — 10 degrees 
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8.3 Ordering Scheme 
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Throughout the formulation of the equations of motion as described above, an ordering 
scheme was applied to reduce the number of terms in the equations. The orders of the various 
variables are as shown here. The largest terms in the equations of motion were of order e. It was 
decided to retain all terms of orders e , e*, e 2 , and e 2 and only linear terms of order e 3 . This 
convention is applied throughout the derivation of the blade and hub matrices. Note that terms 


such as <5/3 refer to the perturbational quantities rather than the variational quantities. 

Order 1 cos cos# 0 i 7> I £ > 

A, /x, i/( 3 , S^, sin sin 6 > 0 , 5 X , , 0 O , # 0 , x x 

tan 63 , t an < 54 , tan < 5 5 

A A A % 49, % C.C.C, «C, *C, <?>C, Sx F , 8y F , Sz Fl 

fig, 8(j > s , a s , 8a s ,ip s , 8'4 > s ,(3qC’> {3gs-> 

■ .. . . / / 

sin cos 0 , 0 , 0 , 0 , 0 , 

V 6 , Wfr, life, v' b , V b , V b , w b , w b , w b , w’ b ,w b , w b , 

y i,2/i,2i,zi 

ftp' jt ’ ~a >Ci r Ir > ^CS> Ccsi i ^-'m. 

*^57 0si ^si ®si &Si ^si .s 1 V’s) $^Psi Ws, 

x F ,x F , 5x f , 6 x f , y F ,y F , 6y F , Sy F , z F , z F ,Sz F , 6z F , 


Order e 


Order e ? 


l,U,U,U e ,U e ,U e ,±i, 

Order e 2 JJ A v b v' b drj d( , f f A v b v b dp d(, [f A v b v b drj d(, f f A v b v b d-q d(. 

If A w b w b dT > d C,.fl A w b ll ’b d n d C, . 1.1 a w'bw'b dr > d (> J.f A drj dQ 



8.4 Coordinate Systems 
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Inertial 


Wing 


Pylon 


Rotating 


Undeformed 


Cross- 

section 


Deformed 


centered 
at wing ea 
(undef) 

centered 
at wing ea 
(undef) 

centered 
at wing ea 
(deformed) 

centered 


pylon CG 
centered 


centered 


centered 


centered 
at x along 
blade 

centered 
at x along 
blade 


_ J j Ground 


j^Ground 

1 ( L 


= P a] 


= [T kl 


= [Trh] 


= {T cu } 


= [T dc 


translation by 
steady 

vehicle motion 

translation by 

( 0 , 0 , 0 ) 

translation by 

(*/i i Vh > %h ) 

offset by 

( x pcg j Vpcg ) Zpcg ) 

offset by 

(xcGiVCG, h) 

translation by 

( 0 , 0 , 0 ) 

translation by 

( 0 , 0 , 0 ) 

translation by 

(M,0) 

translation by 

(■U, v c , w c ) 


rotation by 

( 0 , 0 , 0 ) ' 

rotation by 

(0,0, A) 


rotation by 

(0,Q p ,0) 

rotation by 

( 0 , 0 , 0 ) ' 

rotation by 

(0,0, VO 


rotation by 

(#o, 0, 0) 

rotation by 


Note that the Body coordinate system as described above is also an inertial frame of 

reference. The transformation matrices as defined above are given as follows: 

cos A sin A 0 
[7a] = -sin A cos A 0 


[T wi ] = 


cos xp s sin xp s 0 

— sin ip s cos xp s 0 

0 0 1 


— sin A cos A 0 
0 0 1 

cosa 5 0 sina 5 
0 1 0 

— sin a s 0 cos a s 


1 0 0 

0 cos <p s — sin <p s 

0 sin <p s cos cf) s 


[Tw] = 0 


cos a p 0 — sin a p 


1 0 

0 cos a. 


M = 


cos xp sin ^ 0 

— sin xp cos xp 0 
0 0 1 
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’ cos P 0 

sin p p 

[Tur] = 

0 

1 

0 


— 

sin P p 0 

cos P p _ 


‘ 1 

0 

0 

[Tcu] = 

0 

COS 00 

sin 6 o 


0 

— sin 0o 

COS 0Q 


1 — 


1 '2 


i w '2 

2 V b 


w' b sin 9q + v' b cos 6$ —v b sin 9q + w b cos 0 O 
cos (p sin <p 

cos <p 


— sirup 


2^6 

, ( — w' b cos cp + v b sin <p) sin 0o 

[^dc] =| — (v b cos <p + iu' b sin <p) cos 0o 

(v' b cos <p + w b sin <p ) sin 9q 
+ (—w' b cos <P + v' b sin <p ) cos 0o 

With the gimbal motion included, the matrix [T ur ] as given above is modified to be 

cos 0 p + sin (3 p (0 GG cos ip + 0 GS sin ip) 0 - cos p p (0 GC cos ip + p GS sin ip) + sin j3 

[Tur] =0 10 

. cos Pp ( Pgc cos + Pgs sin VO - sin (3 p 0 cos (3 p + sin p p {(3 GC cos ip + fi GS sin ip) 
With the couplings included, (i.e. 9 = 9- k p/3 p - k ^ + k p p G (0 GC cos ip + /3 GS sin ip)) 


these matrices change as shown below: 


where 


1 0 

0 ( 1 + T cu ) cos + T'L sin 0 O 

0 T^cos0 o - (l + 7l)sin0 o 


0 

~T^ U cos 0o + ( 1 + T^) sin 9q 
(l + T, cos 9 0 + T\ ^ sin 0 O 


t L — —kpff ft kpc. C + ( 0 4 - kp^ Q ( k p p G 8 GC cos ip + ^p8G Pcs sin ip) 
Tcu = kpP P + kpC C — k v p G p GC cos ip — k p p G p GS sin ip 


+k p 0 0 k p ( C (k p 0 G 8 ac cos 0 + k p p G 8 as sin ip) 


rp\2 
1 dc 



w b sin Oq 

-v b sin 0Q 

_ L w ’i _ i ? /2 

2 w b 2 t> 

+v b cos #o 

+w b cos do 
u T h 
^ 1 dc 

(— w b cos <p + v b sin (j>) sin 9q 

+^c 2 

— (v b cos <p + w b sin (j) ) cos #o 
+Tdc 

(v r b cos (p + w' b sin <f>) sin 0q 

cos<p 

sin 4> 

+ (— iv b cos (j> + v b sin <fi) cos Oq 

+T£ 

— sin 0 

COS (j) 



( kpQ kpflG kp(3 kp£ kpj3G ) W b 

Skp^kppQ — kpp — kpgQ — kp^ — Skppkp^ + 3 kppkppQ 

kppkp<;kppG + 2 kppkp^kppG ~ ~kppk p Qkpf3G ~ -kppk^kppc 


cos 6q 
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+ 


( kp^ k p 0G 4~ kpp kpftG kppkp^) w b 
~ kpfiG 4~ kppkpC'kpfiG 4~ krf-kppc + ^pC — k pP k P< 

+ I ~~ k p£ k pfi ~4 k p fi 4 ~ ^kpfikp^kpflG “ k pC k pj3G 

kpfi ~ k pfi k p£ k pfiG 4" ^p/3^^p/3G 4~ k pfik p fiG 

kpfi ~ k p8G — k lc ~ ^ k p& k p(, + 3 k p(3 k pj3G 


-k 2 
k vfiG 


sin^o 


3kry{ kpQG 

r> rr'-' jy^L/V-r pv, /~r~ /“ "* f r- j-j 

~ k p3 k p^ k ppG ^ k pP k pC k p(3G “ 2k 2 p k p^ k p pG — Zkppkp^kpfiG 


• 1 « J 4 ' 

' (^p£ “I” ^’p /3 4 “ k pf3^p(, k p8G kpflG) 


w\ 


COS 00 


+ 


V/3G 


- K p(3G 


k^k 


'Pi^pSG kpC^pfiG 


j pfi 4~ k p(3 4~ 4 kpfi k p<^ k p/3G ^pC 


+ ^p/?^pC 

f A^n 


^P /3 k 2 fik^kpfiQ ^p / 3 


k 2 

- k pPG 

k k k lpG 4 * k p{3 k PPG 


Wu 


4- 


+ 


+ 


rr - • r > ri rr 

4 “ ( A^p /3 A-'p<; kppkppc A’p^ A^p^(^) V b 

(kp$ + fcp^ 4- kppkp^kpflG} w b 

k p{3G 4“ ^p/3 4~ ^kppk^-kppG + 3/Cp^/cp^ — 3 kp^kppc - 

+2k 2 p kp( kpp Q + ^p^^^p/ 3 G — 3 k pfi k p{3G + ^ 

( k p fi k pC, kpfikpfiG A^p^/Cp^c?) 

+ ^p^G^PC ~ k Pfi k pQ k p(3G + ^p/3^p < k pfiG 
k pj3 k p3G + A:^fc pC 4- k^kpp 


sin 0o 


2 k pfi k p{ k ppG 


COS 00 cos 0 


4 A^£? Ajp^ kpfiG 


+kl 


: pfiG k pfi 


k p/3 k pC k pj3G 


k p^ k pfiG kpfi 4 " kppG A* p £ 


sin 0o cos (f) 


kpfiG 4~ ^p/3 4“ Zkppk^kppG + 3 A Ipfikp^ ^kp^kppG ~ k pP k pC k pj3G 

J r2k 2 pk p £kp{3G 4- k pp k ^ k ppG “ ^ k p/3 k pPG 4~ k 2 ^ 


4 - ( kpp 4 " k p pQ kppkp^kpQG k pk) ^ b 


+ 


4 kpp kp^ k p p g 4* k 2 


p/3G k pC k PP k p<^ k pf3G k ^p k pQ k PP G 2 + k p8G k 
k p3 k p( k pi3G k p0 k pPG + kppkpQ -4 kp^ k p p — k^kpfiG 

kpfi 4~ kpfiG k pC, 

4“ ( kp^ kpfiG 4~ kpfi kpfiG kpfikp v b 


>pfi 


Wu 


Wu 


cos 0 q sin <f) 


sin 0o sin cf> 


k 2 

k pPG 


k 2 

K P fi 


2 kpfi k 2 s kpfiG 4~ 3 kpfikp^ 3A ^ kpfiG -kpfikp^k pG 

> ■, . s , 0 , 0,0 ^ , , , 2 ‘ 

; pC 


+ < 2kpfik p ( t k p fiG 


+ kppk^-kppQ 3 k 


'ppkppG + k? 


w. 


kppkrfkppG - k p v' b 


+ 


+ 


+ 


4" {—kpfi 4- kpfiG 
—Akpfikp^kpfiG 4- kpfigkp^ — kpfik 2 ^k 2 fi G 4- A^A^Arp^c 
+ kpfiG k pfi — kpfikpC k pfiG ~ k pfi k pfiG 4" ^p^^pC 4" k p^ k p fi 
-k^kpfiG - A: p ^ + k p fiG 


cos 0o cos <p 


Wu 


k 


pC 


sin 0 q cos <fi 


4- 


+ 


4“ {kp^ kpfiG 4~ kpfikpfiG kpfi k p ^ ) v b 

( kpfi 4 " kpfiG kpfikp^ k p fiG A^) 

—3kppkpQ — A;^ — A:^^ k 2 ^ k 2 ^ G — 2k p fik^k p fiG — 2k 2 fikp^k p fiG 
+2k p fikp£k p p G — k p p G — k 2 fi 4- 3k p fik p fiG 4- 3kp^k p fiG 


cos 0 q sin (f) 


%fi 


{kp^kpfiG 4~ kpfikpfiG kpfikp ^ ) t/a 
k r)3 k pc 4“ ^kpfikpf kpfiG 4~ kpr -4 k ^kp 


^p/3^PC 4“ ^^pfitZp^KpfiG ~t~ ^p^ “4 fcp^kpfiG 
-kpfiG k pP 4- k 2 fikpQk 2 p G — k^kpfi + k 2 fik p an — A;? 

4 “^p^^pC kpfiG “ kpfik^kpfiG — A:. 


^p/?G 

^p/3G 


'pfiG k pQ 


sin 0o sin^ 



8,5 Rigid Blade Simplification 
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The following substitutions are made: 

u e — 0=0 

Vb = —sin £ (x — e*) 

= sin/5(x — ey) 

Therefore, the following derivatives with respect to time and length out the blade are also 

zero: 

» M . / f 

U e U e 0 0 0 0 0 

The non-zero derivatives are: 

Vb = — C cosC (x — ei) wt, = (3 (x — ej) cos/5 

Vb = C ( x - e f) sin C ^ C i x - e i) c °sC = -/? 2 (x — e/) sin/5 + /5(x — e/) cos/5 

v b = — sin C = sin /5 

= ~C cos C w b = 0 cos /5 

■*' • 2 / . 2 

i>b — C sin C — C cos C = — j3 sin (3 + (3 cos (3 

The cross-sectional integrals arising from the foreshortening terms in the kinetic eneigy 
are simplified as follows: 

II A v b v b d V dC, = x sin 2 C, f f A w b w b drj d( = x sin 2 (3 

n A v h v b drj d £ = x£ sin £ cos £ j f A w b w b dr] = x(3 sin (3 cos /5 

J/a v'bV b d v = - c C cos 2 C [J A w b w b dr] d( = xff cos 2 6 

.1.1 A v b v b (lr I d C = sin C (C cos C “ C 2 sin 0 Jj A w’ b w' b dr] d(; = sin (3 (ft cos j3 - if sin 0\ 

It is important to note that all differentiation is completed before any small angle 

assumptions are made. Once these substitutions are made, though, /5 and £ are assumed to be small 



so that 
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v b = -< (x - et) w b = 0 (x - ej) 

v b = ~C ( x ~ w b = 0( x - e f ) 

Vb = C 2 (x - e/)C - c(x - e<) ri)(, = -<? 2 (x -e f )0 + 0(x- e f ) 

v b = ~C = 0 

v b = ~C = (i 

^ = c 2 c - c «>; = -/3 2 /3+/? 

If A v b v b dr l d C = x ( 2 ff A w'bw'b dr ) d( = X/3 2 

If A V b i ’b dr l d C = X (C ff A w b w b dr/ d( = x00 

If A KK dr l d ( = X C If A w'bw'b dr) d< ^ = X @ 2 

If A v bK dr l d C = C (c ~ C 2 c) If A w' b w b drjdC = 0(0- 0 2 0^ 



8.6 Wing Element Matrices 


79 


From the kinetic energy formulation for the wing, a structural mass matrix is calculated. 


[M 5 ] = 


[M?i] = 


[A4] = 


[Mir] = 


[• A4 ] = 


M n 

g AL 

3 

0 

0 

0 

0 

0 

pAL 

6 

0 

0 

0 

0 

0 

g AL 

6 

0 

0 

0 

0 

0 

g AL 

3 

0 

0 

0 

0 

0 


m? 2 ■ 

_ 

0 

UpAL 
35 , 

11 gALA 
210 
0 

0 

0 

0 

OpAL 
70 , 

13 pAL 2 
420 
0 

0 

0 

0 

9pAL 

70 

\3pAL 2 

420 

0 


0 

UpAL 2 

210,, 

g A±A_ 

105 

0 

0 

0 

0 

13 pAL 2 
420, 
pALA 
140 
0 

0 

0 

0 

UpAL 2 

420 

pAL' 

140 

0 


0 

0 

0 

ISpAL 
35 , 

11 pAL 2 
210 

7pALy ( 'M 

20 

o 
o 
0 

9gAL 

1 sJ/Il* 


0 

0 

0 

\3pAL 

35 , 

UpAL 2 
210 
0 

0 

0 


0 
0 

0 

11 pAL 1 

pB 

105 

0 

0 

0 


420 
.\pALijrM 
20 

0 

0 

0 

9pAL 

70 

!3pAL~ 

420 

3pALy C 'M 


o 
0 
0 

UpAL' 2 

Ml 

105 

pAL 2 yc^, 

20 

0 

0 

0 

UpAL 2 


20 


0 

0 

0 

7p/lLvr-A> 

.. 2,0 


0 

0 

0 

13ML 
35 , 

ilpAL 2 

210 

lpALy,-M 

20 


140 
gA LAjjrj, 

30 

0 

0 

0 

UpAL 2 
420 , 
pAL' 
140 

pAL 2 yc\t 

30 

0 

0 

0 

UpAL 2 

P B 


o 

o 

0 

3pALyc'\t 

20 

pAL Vc ’ a / 


105 

pAL 2 yc\t 

20 


6 

0 
0 
0 

3pALyr\t 
20 

pAL 2 y ( 'M 

pM 

6 

0 

0 

0 

7pALy ( \\i 
20 

pAL 2 y r .\t 




The strain energy formulation gives the structural stiffness matrix. 


[K 8 ] 




K?2 ' 





EA 

L 

k! 2 j 
0 

0 

0 

0 

0 

0 

11 ML. 
qmL. 

O ia 

qKL 

X 2 

0 

0 

0 

0 

a ML 

4 L 

0 

0 

0 

0 

0 

0 

u 5a 

-EH,, 

6 ~v 


0 

0 

0 

0 

4# 

0 

0 

0 

0 

0 

0 

GJ 

L 
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[*fa] = 


V/' kl 


[*&] = 


0 

- 12 % 


0 

~q£L. 

V T 2 


The couplings in the wing are added into the stiffness matrix. 


{K c } = 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A46-L 

0 

0 

0 

0 

0 

-A 46 i 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A 

0 

0 

0 

0 

0 

-A «i 

0 

0 

luez 

0 

A45J- 

0 

0 

0 

— A46X 

0 

-A45* 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A467J 

0 

0 

0 

0 

0 

A 46^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

~A 4 5X 

0 

0 

0 

0 

0 

A45 i 

0 

0 

-A 

0 

-A45 i 

0 

0 

0 

A 46 i 

0 

k 45 i 

0 


The work energy from the aerodynamics gives the motion-independent forcing vector and the 



motion-dependent stiffness and damping matrices. 


= 

w.i = 

m = 

[^i = 
[Ki y = 


{F*} = 


-pV z cc L 


1 2 

7>P V CC L 

0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 

0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 

0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 

0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


k a 


1 1 


r 


a 12 

k a 
A 22 

0 
0 
0 

^ sin A 


.2 

sin AL 


jx sm 
^ sin Ae 

0 

0 

0 

— ^ sin A 

— jrj sin A L 

— | sin Ae 

0 

0 

0 

4 sin A 
tq sin AL 
| sin Ae 

0 

0 

0 

— sin A 
sin AL 
— -jsinAe - 


0 

0 

0 

7« 0 L 

jW 2 

^ea 0 L 

0 

0 

0 

^OCqLt 

— J2 a 0^ 

^eaoL 


0 

0 

0 

- sin AL 

0 

-j 2 sin A Le 

0 

0 

0 

sin AL 
A sin AL 2 
^ sin A Le 

0 

0 

0 

-A, sin AL 

— sin AL 2 
•pj sin A Le 

0 

0 

0 

-jq sin AL 
0 

~ sin A Le 


0 

0 

0 

^ cos AL 
^ cos AL 2 
^ecos AL 

0 

0 

0 

cos AL 
^ cos AL 2 
|ecos AL 


2(j) * 
TiO 


_7_ 

2 ? 


' 12 
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2 P v CC L 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13 L 
35 V 

11 L 1 
210 V 

0 

0 

0 

0 

9 L 
70 V 

0 

0 

0 

11 L 2 

1 L* 

0 

0 

0 

0 

13 L 2 

210 V 

105 V 

420 V 

0 

0 

0 

7 eL 
20 V 

1 eL 2 
20 V 

0 

0 

0 

0 

3 eL 
20 V 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 L 
70 V 

13 L 2 
420 V 

0 

0 

0 

0 

13 L 
35 V 

0 

0 

0 

13 L 2 
420 V 

1 L X 
140 V 

0 

0 

0 

0 

11 L 2 
210 V 

0 

0 

0 

3 eL 
20 V 

1 eL 2 
30 V 

0 

0 

0 

0 

7 eL 
20 V 


0 

0 

0 

13 L? 
420 V 
1 L :i 
140 V, 

I eL J 
30 V 

0 

0 

0 

II L? 

210 V , 

1 _ Li 

105 V 
1 eLj 
20 V 


0 

0 

0 

0 

0 


0 

0 

0 


O O O O 



8.7 Assembly of Wing and Pylon 
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The [A/], [C], and [K\ from the coupled rotor hub equations of motion are transformed into 
the fixed frame using a Fourier Coordinate Transformation. Tis transformation takes the degree of 
freedom vector 

[ 00 01c 01s Co Clc Cls ••• X F VF Z F a s <t> s Xp s f3 GC 0 GS J 

and returns a new degree of freedom vector of 

|_ 00 0\c 01 s *•* Co Clc Cla ••• u V N~ hi V N+ 1 W N+ 1 Wjv+1 ^Af-f-2 0GC 0GS J 

Then these matrices are constrained to remove the u degree of freedom. Note that this is only true 
when modeling a wing rigid in extension. This constraint will also be removed when fuselage 
motion is included. After this constraint, these matrices can be added into the wing tip degrees of 
freedom. 

[ 00 01 c 01 s - - * Co Clc Cl a V N+ 1 tJ AT+l ^JV+1 w jy+i <Aat+2 0GC 0GS J 

Note that the (f> degree of freedom used here is <frw +2 or the twist on the pylon side of the downstop 

spring. The wing tip will be 0, the twist on the wing side of the downstop spring. 
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[Tpostj 


T 5 


ZQFF 




1 


0 


0 




0 

0 


0 



0 

0 






0 


1 


0 




0 

0 


0 



0 

0 






0 


0 


1 




0 

0 


0 



0 

0 





_ 






2N h x2N, 





J 

2Nb x6 



2 N h 

x2 



0 


0 


0 ; 










: 0 

0 = 






0 


0 


0 










0 

0 













T 5 













0 


0 


0 










0 

0 





_ 






6x2N t 











6x2 




0 


0 


o:i 




0 

0 

. . . 

0 ' 



1 

0 1 






0 


0 


0 




0 

0 


0 



0 

1 




*- 
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voff = vcg + ypcg 
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The degrees of freedom from the rotor / hub equations of motion include hub motion 
already written in the inertial frame, i.e. after the shaft tilt due to pylon rotation. However, the 
pylon rotation must also be taken into account for the pylon masses and inertias. Beyond this 
difference, the pylon inertias are also centered at the pylon center of gravity rather than at the hub. 
This gives different values for the offsets as shown below. 
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